Consistency of LEP Z + 6-jets excess with an /i — ► aa Decay Scenario 
and Low-Fine- Tuning NMSSM Models 



Radovan Derrmsek* and John F. Gunion^ 

* School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 
^Department of Physics, University of California at Davis, Davis, CA 95616 



in 
o 
o 

(N 

O 

<D 

Q 

\o 

> 

(N 
(N 

m 
O 

in 
o 



i 

a: 

CD . 

43 ■ 



X 



We examine the LEP limits for the Zh — > Z + 6's final state and find that the excess of observed 
events for mj, ~ 100 GeV correlates well with there being an m;, ~ 100 GeV Higgs boson with SM- 
like ZZh coupling that decays partly via h — > bb+r + T~ [with B(h — > bb) ~ 0.08] but dominantly via 
ft — > aa [with B(/i — > aa) ~ 0.9], where m a < 2ms so that a — > t + t~ (or light quarks and gluons) 
decays are dominant. Scenarios of precisely this type arise in the Next-to-Minimal Supersymmetric 
Model for parameter choices yielding the lowest possible fine-tuning. 



LEP has placed strong constraints on Zh production 
where h decays primarily to 6-quarks. The limits on 
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[9zzh/9zzh s J B ( h bb ) arc shown in Fi g-IIl El- 



These limits rule out a Standard Model (SM) Higgs boson 
with mh ^ 114 GeV at 95% CL. However, the plot also 
exhibits the well-known excess of observed vs. expected 
C|jj limits for a test Higgs mass of mh ~ 100 GeV. 
This excess is particularly apparent in the 1 — CLj, result 
(Fig. 7 of p|) obtained after combining all four LEP ex- 
periments. Various interpretations of this excess in terms 
of a non-SM Higgs sector have been suggested 0,0- ^ n 
this letter, we point out that this excess is consistent with 
a scenario in which the Higgs boson has SM-like ZZh 
coupling, but has reduced B(h —> bb) by virtue of the 
presence of h decays to a pair of lighter Higgs bosons, 
h — > aa, where we use the notation a appropriate to 
the NMSSM cases of interest where a is a CP-odd Higgs 
boson. (For a generic two-Higgs doublet model with or 
without extra singlets, the light Higgs could alternatively 
be CP-even, or if CP-violation is present, of mixed CP- 
nature.) As we shall discuss, the a must not have signifi- 
cant B(a — > bb). The most naturally consistent possibil- 
ity is m a < 2rrib so that a — > t + t~ or light quarks and 
gluons. We then emphasize that parameter choices for 
the Next to Minimal Supersymmetric Model (NMSSM) 
that yield the smallest possible fine-tuning typically pre- 
dict precisely this kind of scenario. 

The limits on C^j j apply when the only decays of the 



h to 6-quarks are direct, h — * 
h to 6's were via Zh — > Zaa 



bb. If the only decays of 
— > Zbbbb, alternative con- 
straints 1 on Cf ff = [g 2 zzh / g 2 zzhsM ]B(h - Zaa)[B[a - 
bb)] 2 would apply. Again, C^// is above the background 
Monte Carlo predictions, in particular for m/, Si 105 GeV 
and m a ~ 30 45 GeV. However, since in the LEP an- 
alyzes for the Z2b and Z\b final states the selected can- 
didates are correlated, the Z2b and ZAb excesses cannot 
be treated as being independent if the Higgs boson de- 
cays to both kinds of final state. In particular, a model 
for which the predictions for and are indi- 

vidually allowed but both are close to their respective 
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: Expected and observed 95% CL limits on C^ff = 
izzh/9zzh SM ]B(h -» 66) from Ref. Q are shown vs. 
rrih- Also plotted are the predictions for NMSSM pa- 
rameter choices in our fixed tan/? = 10, Mi t 2,3(mz) = 
100, 200, 300 GeV scan that give fine-tuning measure F < 25 
and m ai < 2n%b and that are consistent with the preliminary 
LHWG analysis code. 

observed exclusion limits is actually ruled out because 
the combined Z2b and ZAb rates give too large a net 
Z + 6's event rate The escape we propose is that 
m a < 2mb so that a — ► t + t~ (or a — > light quarks 
and gluons if m a < 2m T ) is dominant, implying that 
Zh — ► Zaa does not lead to a Z + 6's final state. If the 
ZZh coupling is full SM strength, then m h ~ 100 GeV 
with B(h -> bb) - 0.08 and B(h -> aa) - 0.9 fits the 
observed Z2b excess nicely. Meanwhile, there are no cur- 
rent limits on the Zh — > Zaa — > Zt + t~t + t~ final state 
for m h > 87 GeV @. Further, since the ZZh and WWh 
couplings are very SM-like such a model gives excellent 
agreement with precision electroweak data. 

We are particularly led to the above interpretation 
of LEP data as a result of our successful search [{| for 
NMSSM parameter choices for which fine-tuning of the 
model with respect to unification scale parameters is 
absent. We found that the NMSSM parameters that 
yield small fine-tuning typically predict an hi and a\ 
with characteristics similar to and sometimes exactly 
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matching the requirements described above. In contrast, 
in the CP-conserving Minimal Supcrsymmctric Model 
(MSSM), h — > bb decays are dominant and all param- 
eter choices consistent with LEP limits on are such 
that the fine-tuning and hierarchy problems are severe. 

The NMSSM is an extremely attractive model [(|. 
First, it provides a very elegant solution to the \i problem 
of the MSSM via the introduction of a singlet superfield 
S. For the simplest possible scale invariant form of the 
superpotential, the scalar component of S naturally ac- 
quires a vacuum expectation value of the order of the 
SUSY breaking scale, giving rise to a value of fj, of order 
the electroweak scale. The NMSSM is the simplest su- 
persymmetric extension of the standard model in which 
the electroweak scale originates from the SUSY break- 
ing scale only. Hence, the NMSSM deserves very serious 
consideration. 

Apart from the usual quark and lepton Yukawa cou- 
plings, the scale invariant superpotential of the NMSSM 
is W = XSH u Hd + fS 3 depending on two dimension- 
less couplings A, k beyond the MSSM. [Hatted (un- 
hatted) capital letters denote supcrficlds (scalar super- 
field components).] The associated trilinear soft terms 
arc \A\SH u Hd + ^A K S 3 . The final two input param- 
eters are tan (3 = h u /hd and /i e ff = As , where h u 

ll : . hd = (Hd) and s = (S). The Higgs sector of 
the NMSSM is thus described by the six parameters 
A , k , A\ , A K , tan {3 , /x e ff ■ In addition, values must be 
input for the gaugino masses and for the soft terms re- 
lated to the (third generation) squarks and sleptons that 
contribute to the radiative corrections in the Higgs sector 
and to the Higgs decay widths. 

The particle content of the NMSSM differs from the 
MSSM by the addition of one CP-even and one CP-odd 
state in the neutral Higgs sector (assuming CP conser- 
vation), and one additional neutralino. The result is 
three CP-even Higgs bosons (^1,2.3) two CP-odd Higgs 
bosons (01,2) and a total of five neutralinos X? 2 3 4 5- The 
NMHDECAY program 0, which includes most LEP con- 
straints, allows easy exploration of Higgs phenomenology 
in the NMSSM. 

In H, we found that the Next to Minimal Supersym- 
metric Model (NMSSM) can avoid the fine-tuning and 
hierarchy problems of the MSSM for parameter choices 
that were consistent with all LEP constraints available 
at the time. Defining the fine-tuning measure to be 



F 



MaxpFp = Maxp 



dlogmz 



dlogp 



(1) 



where the parameters p comprise all GUT-scale soft- 
SUSY-breaking parameters, we found that F < 10 could 
be achieved. Further, at moderate tan/3, the parame- 
ter choices with such F are always such that the light- 
est CP-even Higgs boson, hi, is SM-like as regards its 
gauge and fermionic couplings, but decays primarily into 
a pair of the lightest CP-odd Higgs bosons of the model, 



hi — > aiai. The importance of such decays was first 
emphasized in |p , and later in H|, followed by extensive 

work in munia. 

We note that a light ai is natural in the NMSSM in the 
kA k , \A\ — » limit. This can be understood as a conse- 
quence of a global U(1)r symmetry of the scalar poten- 
tial (in the limit kA k ,XA\ — > 0) which is spontaneously 
broken by the vevs, resulting in a Nambu-Goldstonc 
boson in the spectrum 0. This symmetry is explic- 
itly broken by the trilinear soft terms so that for small 
kA k ,\A\ the lightest CP odd Higgs boson is naturally 
much lighter than other Higgs bosons. For the F < 10 
scenarios, \(niz) ~ 0.15 -j- 0.25, n(mz) ~ 0.15 ~ 0.3, 
\A K (m z )\ < 4 GcV and \A x (m z )\ < 200 GeV, imply- 
ing small kA k and moderate \A\. The effect of \A\ 
on m ai is further suppressed when the ai is largely sin- 
glet in nature, as is the case for low-F scenarios. There- 
fore, we always obtain small m ai . We note that small 
soft SUSY-brcaking trilinear couplings at the unifica- 
tion scale are generic in SUSY breaking scenarios where 
SUSY breaking is mediated by the gauge sector, as, for 
instance, in gauge or gaugino mediation. Although the 
value A\(mz) might be sizable due to contributions from 
gaugino masses after renormalization group running be- 
tween the unification scale and the weak scale, A K re- 
ceives only a small correction from the running (such 
corrections being one loop suppressed compared to those 
for Ax). Altogether, a light, singlet ai is very natural in 
models with small soft SUSY-brcaking trilinear couplings 
at the unification scale. Finally, we note that the above 
\(mz) values are such that A will remain perturbative 
when evolved up to the unification scale, implying that 
the resulting unification-scale A values are natural in the 
context of model structures that might yield the NMSSM 
as an effective theory below the unification scale. 

We will now discuss in more detail results for the 
NMSSM using the representative fixed values of tan /3 = 
10 and Mi l2 , 3 (m z ) = 100,200,300 GcV while varying 
all other model parameters. Similar results are obtained 
for other choices of tan/? and M\2z(rrhz)- The points 
plotted for the NMSSM in Fig. Q] show the Cf ff pre- 
dictions for all parameter choices in our scan that had 



F < 25 and 



< 2mt> and that are consistent with 



the experimental and theoretical constraints built into 
NMHDECAY as well as with limits from the prelimi- 
nary LHWG full analysis code. The eight F < 10 points 
are singled out. From Fig. ^ we see that these latter 
points cluster near m/^ ~ 98 ~ 105 GcV (see also Fig. 3 
of 0). We will sec that most are such that and 
B{h\ — ► bb) are appropriate for explaining the Cfff ex- 
cess. The other primary h\ decay mode for all the plotted 
points is hi — > ai<zi with ai — ► t + t~ or light quarks and 
gluons (when m ai < 2m T ). In TableHJ we give the precise 
masses and branching ratios of the hi and ai for all the 
F < 10 points. We also give the number of standard devi- 
ations, n b s («oxp) 7 by which the observed rate (expected 
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TABLE I: Some properties of the hi and ai for the eight 
allowed points with F < 10 and m ai < 2mb from our 
tan/3 = 10, Mi, 2 , 3 (mz) = 100,200,300 GeV NMSSM scan. 
The n bs, n oxp and s95 values are obtained after full process- 
ing of all Zh final states using the preliminary LHWG analysis 
code (thanks to P. Bechtle). See text for details. Ng^ c is 
the statistical significance of the best "standard" LHC Higgs 
detection channel for integrated luminosity of L = 300 fb _1 . 

rate obtained for the predicted signal+background) ex- 
ceeds the predicted background. These are derived from 
(1 - CL^obscrvcd and (1 - CL 6 ) oxpGCtod using the usual 
tables: e.g. (1 - CL b ) = 0.32, 0.045, 0.0027 correspond 
to lcr, 2a, 3a excesses, respectively. The quantity s95 is 
the factor by which the signal predicted in a given case 
would have to be multiplied in order to exceed the 95% 
CL. All these quantities are obtained by processing each 
scenario through the full preliminary LHWG confidence 
level/likelihood analysis. If n cxp is larger than n b s then 
the excess predicted by the signal plus background Monte 
Carlo is larger than the excess actually observed and vice 
versa. The points with < 100 GeV have the largest 
n b s ■ Point 2 gives the best consistency between n Q b s and 
^oxpj with a predicted excess only slightly smaller than 
that observed. Points 1 and 3 also show substantial con- 
sistency. For the 4th and 7th points, the predicted excess 
is only modestly larger (roughly within la) compared to 
that observed. The 5th and 6th points are very close to 
the 95% CL borderline and have a predicted signal that is 
significantly larger than the excess observed. LEP is not 
very sensitive to point 8. Thus, a significant fraction of 
the F < 10 points are very consistent with the observed 
event excess. 

We wish to emphasize that in our scan there are many, 
many points that satisfy all constraints and have m ai < 
2m;,. The remarkable result is that those with F < 10 
have a substantial probability that they predict the Higgs 
boson properties that would imply a LEP Zh — > Z + 6's 
excess of the sort seen. The smaller number of F < 10 
points with m ai substantially above 2rrib all predict a net 
Z + 6 ! s signal that is ruled out at better than 99% CL by 
LEP data. Indeed, all F < 25 points have a net h — > 6's 
branching ratio, B(h\ — > 66) + B(h\ — > aidi — ► 6666) > 
0.85, which is too large for LEP consistency if m ai is 
substantial. Analysis of points with m ai very near bb 
decay threshold, but such that a\ — > 66 is dominant, is 



very subtle. Such points arise for F < 10 and require 
further analysis in cooperation with the LHWG. 

An important question is the extent to which the type 
of h — > aa Higgs scenario (whether NMSSM or other) 
described here can be explored at the Tevatron, the LHC 
and a future e + e _ linear collider. This has been exam- 
ined in the case of the NMSSM in @, E3, EH , witn thc 
conclusion that observation of any of the NMSSM Higgs 
bosons may be difficult at hadron colliders. At a naive 
level, the hi —> a\ai decay mode renders inadequate the 
usual Higgs search modes that might allow hi discov- 
ery at the LHC. Since the other NMSSM Higgs bosons 
arc rather heavy and have couplings to 6 quarks that 
are not greatly enhanced, they too cannot be detected at 
the LHC. The last column of Table |U shows the statisti- 
cal significance of the most significant signal for any of 
the NMSSM Higgs bosons in the "standard" SM/MSSM 
search channels for the eight F < 10 NMSSM parameter 
choices. For the hi and ai, the most important detection 
channels are hi — > 77, Whi+tthi -^jj^X^thi/cn -» 
tibb, bbhi/ai — > bbr + r~ and WW — » h\ — ► t + t~ - see 
[l2j |. Even after L = 300 fb _1 of accumulated luminos- 
ity, the typical maximal signal strength is at best 3.5<r. 
For the eight points of Table [I] this largest signal de- 
rives from the Whi + tthi — > 77£ ± X channel. There is 
a clear need to develop detection modes sensitive to the 
dominant hi — > aiai — > t + t~t + t~ decay channel. 

Let us consider the possibilities. One detection mode 
that can be considered is WW — > hi — > a±ai — > At. 
Second, recall that the X2 ~ ¥ hiXi channel provides a 
signal in the MSSM when hi — ► 66 decays are domi- 
nant. See, for example, ^3|- It nas n °t been studied 
for hi — * a\di — > 4r decays. If a light Xi provides 
the dark matter of the universe (as possible because of 
the y? y? — > ai — ► X annihilation channels for a light 
01, see |l4l IT5I and references therein), the m~ — m~ 
mass difference might be large enough to allow such de- 
cays. Diffractive production [ltH ]. pp —> pphi — > ppX, 
where the mass Mx can be reconstructed with roughly 
a 1 — 2 GeV resolution, can potentially reveal a Higgs 
peak, independent of the decay of the Higgs. A study 
[17] is underway to see if this discovery mode works for 
the hi — > aidi — > 4r decay mode as well as it appears 
to work for the simpler SM Ksm — * bb case. The main 
issue may be whether events can be triggered despite the 
soft nature of the decay products of the r's present in X 
when hi aidi — » 4t as compared to ft-sM - * 66. 

At the Tevatron it is possible that Zhi and Whi pro- 
duction, with hi — > aiai — ► 4r, will provide the most 
favorable channels. If backgrounds are small, one must 
simply accumulate enough events. However, efficiencies 
for triggering on and isolating the 4r final state will not 
be large. Perhaps one could also consider gg — > hi — > 
aiai — > 4r which would have substantially larger rate. 
Studies are needed. If supcrsymmetry is detected at the 
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Tcvatron, but no Higgs is seen, and if LHC discovery of 
the hi remains uncertain, Tcvatron studies of the 4r fi- 
nal state might be essential. However, rates imply that 
the hi signal could only be seen if Tevatron running is 
extended until L > 10 fb _1 has been accumulated. Of 
course, the LHC would observe numerous supersymmetry 
signals and would confirm that WW — > WW scattering is 
perturbative, implying that something like a light Higgs 
boson must be present, but direct detection of the hi 
might have to rely on an extended Tevatron run. 

Of course, discovery of the hi will be straightforward at 
an e + e~ linear collider via the inclusive Zh — > £ + £~ X re- 
constructed Mx approach (which allows Higgs discovery 
independent of the Higgs decay mode) . Direct detection 
in the Zh — > ZAt mode will also be possible. At a 77 
collider, the 77 — > h — > 4r signal will be easily seen |l8j . 

In contrast, since (as already noted) the ai in these 
\ow-F NMSSM scenarios is fairly singlet in nature, its 
direct (i.e. not in hi decays) detection will be very chal- 
lenging even at the ILC. Further, the low-F points are all 
such that the other Higgs bosons are fairly heavy, typi- 
cally above 400 GeV in mass, and essentially inaccessible 
at both the LHC and all but a > 1 TeV ILC. 

We should note that much of the discussion above re- 
garding Higgs discovery is quite generic. Whether the a 
is truly the NMSSM CP-odd ai or just a lighter Higgs 
boson into which the SM-like h pair-decays, hadron col- 
lider detection of the h in its h — > aa decay mode will 
be very challenging — only an e + e~ linear collider can 
currently guarantee its discovery. 

In conclusion, we reemphasize that the prominent LEP 
event excess in the Z+b's channel for reconstructed Higgs 
mass of nth ~ 100 GeV is consistent with a scenario 
in which the ZZh coupling is SM-like but the h decays 
mainly via h — > aa — > t + t~t + t~ (requiring m a < 2mb) 
leaving an appropriately reduced rate for h — > bb. Wc 
strongly encourage the LEP groups to push the analysis 
of the ZAt channel in the hope of either ruling out the 
h — > aa — > 4t scenario, or finding a small excess consis- 
tent with it. Either a positive or negative result would 
have very important implications for Higgs searches at 
the Tevatron and LHC. Further, we have emphasized that 
the NMSSM models with the smallest fine-tuning typi- 
cally predict precisely the above scenario with h = hi and 
a = a\. We speculate that similar results will emerge in 
other supersymmetric models with a Higgs sector that is 
more complicated than that of the CP-conserving MSSM. 
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